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Abstract 

Recently, fractional derivatives have been employed to analyze various systems in engineering, 
physics, finance and hidrology. For instance, they have been used to investigate anomalous diffusion 
processes which are present in different physical systems like: amorphous semicondutors, polymers, 
composite heterogeneous films and porous media. They have also been used to calculate the heat 
load intensity change in blast furnace walls, to solve problems of control theory and dynamic prob- 
lems of linear and nonlinear hereditary mechanics of solids. In this work, we investigate the scaling 
properties related to the nonlinear fractional diffusion equations and indicate the possibilities to 
the applications of these equations to simulate the water transport in unsaturated soils. Usually, 
the water transport in soils with anomalous diffusion, the dependence of concentration on time 
and distance may be expressed in term of a single variable given by X q = x/t q . In particular, for 
q = 1/2 the systems obey Fick's law and Richards' equation for water transport. We show that a 
generalization of Richards' equation via fractional approach can incorporate the above property. 
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INTRODUCTION 



In these last decades we have witnessed a growing use of the fractional calculus to describe 
diverse systems in engineering, physics and finance. In physics, the fractional diffusion 
equations have been used to describe anomalous diffusion processes which are present, for 
instance, in amorphous semiconductors Q], polymers composite heterogeneous films [3] 
and porous media 4J]. In engineering, for example, the fractional derivative has been used in 
viscoelastic models which generalize the classical mechanical models p. In finance market, 
the Gaussian assumption for the probability distribution of returns has been extended by 
using non-Gaussian assumption and the fractional derivative approach, according to the 
empirical studies Q]. 

In hydrology, the description of the transport of water in many soils, in horizontal column, 
is done by the nonlinear diffusion equation pj 



(1) 
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In this equation presents a scaling dependence of concentration on time t and distance x 
given by 



x 



Ai/a(0) = - • (2) 
This is so called the Boltzmann variable. Eq. written in term of A1/2 is given by 
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However, in the soils like Columbia silt loam, Hesperia silt loam, Hesperia sandy loam, 
Salkum silty clay loam and kaolinite-sand mixtures, deviations from J2I) have been observed 

n 

and the dependence of concentration on t and x |9( is given by 



\(0) = I , (4) 

where q can be greater or less than 1/2. In order to obtain the generalized scaling law 
P|, based on the nonlinear diffusion equation, further assumptions have been incorporated 
into Eq. (0): In 0, the diffusion coefficient is assumed to have an explicit dependence 



2 



on distance given by x~ l3 D{9) ) while in the diffusion coefficient is assumed to have an 
explicit dependence on time given by t 2m ~ 1 D(6), where m is a positive constant. Indeed, the 
above assumptions are phenomenological, and other possibilities to obtain the scaling law 
(J3J) exist. An interesting possibility is to use the fractional approach. Fractional derivatives 
can be defined in terms of integral, and due to this fact they are related to the memory 
effects when they are employed in differential equations applied to physical systems. In 
stochastic processes the fractional derivatives are related to the non-Markovian nature of 
the systems, for instance, they can provide a dynamic framework to describe long random 

n 

walks and waiting time between successive walks 11]. 

In this work, we investigate the scaling dependence on time and distance in nonlinear 
fractional diffusion equations. To do so, we employ the Riemann-Liouville and Caputo defi- 
nitions for fractional derivatives. Indeed, there are other definitions for fractional derivatives 



12j, but the Riemann-Liouville and Caputo definitions are most often used in the literature. 
According to the previous discussions, these nonlinear fractional diffusion equations, with 
the scaling dependence on time and distance, may be useful to describe the water transport 
in soils when the memory effects are present in these systems. The presentation of this 
work is outlined as follows. In section II we discuss on the scaling dependence on time and 
distance in several nonlinear fractional diffusion equations by using the fractional derivatives 
on distance and time. In section III, we point out on some possible applications of these 
equations in the water transport in soils. In section IV, we present our conclusions. 



NONLINEAR FRACTIONAL DIFFUSION EQUATIONS AND SCALING DEPEN- 
DENCE ON TIME AND DISTANCE 

Several physical systems are intrinsically structured with some fundamental scaling re- 
lations. When they are revealed they can bring a beneficial insight into a new analysis of 
the systems. Its presence in physical systems indicates that the quantities involved in the 
scaling relation are more rigidly connected, even with the changing of the parameters of the 
systems. For instance, many systems of the water transport in soils are described by Eq. 

which is related to the dependence of the concentration on distance and time (jBJ)- The 
propagator on fractals in Sierpinski gasket also presents a scaling dependence on distance 



and time in its asymptotic form 



111 ]. A remarkable presence of scaling relation appears in 
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the Gibbs potential on the region of the critical point Q|. 

In order to investigate the scaling relations of nonlinear fractional diffusion equations, we 
consider the Riemann-Liouville and Caputo definitions for fractional derivatives which are 
given by 

1 d n r z fir) 
Dr L(z J(z) = / 7 \ 'i dr (5) 

and 

where D^ L ^ indicates the Riemann-Liouville operator and D^,,, indicates the Caputo op- 
erator. These two definitions are connected by the following expression 

n-l 

D a RL{z) f(z) = D« {z) f(z) + Y,f {k \0 + )<S>(k- a+ i)(z) , (7) 



where 



fe=0 



Z > 



Mz)=\ m ' ' \ ■ (8) 
y 0, z < I 

First, we consider the Riemann-Liouville fractional derivative on distance applied to Eq. 

f = ^x ( ,)[^)£], 0<«<1. (9) 

This equation, with D a {6) = D, has been used to study anomalous diffusion for Levy flights. 
For D a {6) = D9 a , it has been studied in [yj. By using the generalized scaling variable (0J 
we obtain from Eq. (0) the following result 

\ q d6 1 d f x D a (6)d T 6 



q t d\ q r(l - a) dx J (x - r)« dr ' (10) 
Introducing the variable u = r/t q into Eq. (fTUj) yields 

Aft 1 

(11) 
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with q = 1/(1 + a). We see that Eq. (JHJ) can be expressed as a function of generalized scaling 
relation (0J). Note that, for a = 1, we recover the Boltzmann variable. 

For Caputo operator on distance, we also use the above procedure and the result is similar 
to Eq. (HI} that is 
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with q = 1/(1 + a). 

For fractional derivative on time, we consider the Caputo operator acting on the left side 
ofEq. ©, 
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dx 



< a < 1. 



(13) 



This equation has been used to study anomalous diffusion for subdiffusive regimes [n| 
Introducing the variable u = -jS^T into Eq. (|T3*j) and using the variable (HJ) we obtain 
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t a T(l - a) J (X g - u)' 
By setting a = 2q, we have 
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Note that the singularity which appears in the variable u = -^p^T for t — > does not 
invalidate the integration of Eq. (|T3*j) due to the fact that the limits of the corresponding 
integral tend to the same values. Therefore, Eq. ()14)) is valid for t > 0. The scaling relation 
related to Eq. (fT3*|) is also corroborated by the analytical result obtained in An 
observation on the integral of the left side of Eq. ((Hj); It does not represent the usual 
definition of fractional derivative (0) due to the fact that the variable 9 depends on both X q 
and integration variable u. For this reason, we have employed the operator O^,^ instead of 

u c(ty 

For Riemann-Liouville operator on time acting on the left side of Eq. (0), the correspond- 
ing fractional equation related to the scaling variable (@J) is given by 
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where q = a/2. In the case of the Riemann-Liouville operator on time acting on the right 
side of Ecl (JTJ), it can be replaced by the Caputo operator on time acting on the left side 
(JT3J), see [11]. The corresponding fractional equation is 
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WATER TRANSPORT IN SOILS AND NONLINEAR FRACTIONAL DIFFUSION 
EQUATIONS 



As we have previously mentioned that many systems with the transport of water in soils, 
in horizontal column, are simulated by Eq. (1) and scaling relation @. However, the scaling 
relation 



can not provide a good agreement for all the experimental results 
We now consider a generalization of Eq. (0) which has been introduced in 
diffusion coefficient dependent of the time given by 



10j with the 
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where (3 is a real parameter and 9 is the volumetric water content. In particular, Eq. (|18|). 
for DjA9) = D and (3 — — 2, introduced by Batchelor, has been used to simulate turbulent 
lows [15| . For water transport in soils, the boundary conditions are usually chosen as follows 



x > 0, 



9 = 9i, t 



(19) 



x = 0, 9 = O , t > 



(20) 



where 0j is the initial water content along the soil column and 9q is the constant water 
content at x — 0. It has been shown that Eq. (|T8|) is related to scaling relation and 
it has also been used to simulate the single-energy gamma-ray attenuation data for Salkum 
silty clay loam jlfl| . 

Now, we employ an anternative route to analyze Eq. (|TB*j) in comparison with the frac- 
tional equation ((T7j). Integrating Eq. (fTHj) on time we have 



J(Iit) _ 9fe „ )= _^ J r4ts^M dT . (21) 

The constant T(l — (3) was introduced by convenience. Introducing the variable r = t — u 
yields 

* d x [Dp{6)d x 6] 

r(i-/3)7 (t- u y 



9(x, t) - 9(x, 0) = W1 m / \; > du . (22) 



Applying the derivative on time we obtain 

89 _ 1 d f d x [Dp(e)d x 6} A 

at ~ r(i-f3)dt J it - u y u ■ (23) 

It is interesting to note that this last equation looks like Eq. (|17p. but they are not the 
same. In Eq. f!23|) the variable 9 depends on x, t and u, while in Eq. f!17|) the variable 9 
depends only on x and r. Operationally speaking, the denominator of the integrand of Eq. 
()17|) . t — t, runs from t to zero, while the denominator of the integrand of Eq. (|21j) . r, runs 
from zero to t; They are the inverse of each other. Eq. (j!7)) and Eq. (|2"T]) are the same for 
a specific case, i.e., if the numerators of the integrands of these equations are symmetrical 
with relation to the point t/2 for all the instant t. 

An other generalization of Eq. (Q) has been introduced in jt| with the diffusion coefficient 
dependent on the distance 



d9_ _ d_ 

dt dx dx |_~ >J v ' dx 
where J is the flux. In particular, Eq. (J2lj) . for Dp(9) = D and (3 = —4/3, introduced by 
Richardson, has been used to simulate turbulent flows Q|. Eq. @ can be obtained from 
Eq. ()24j) by taking the early contributions of the flux. 

We should mention that the complexities of the soil structures are not considered in the 
descriptions of the water transport in soils by using diffusion equations. In these descriptions 
are considered the changing of porous media and bulk density, and the systems may be 
categorized as follows (1) rigid porous media with no change in the bulk density. (2) 
semirigid porous media with no change in the bulk density and (3) swelling porous media 
with the change in the bulk density. The first category is usually described by Eq. with 
the diffusion coefficient depends only on the volumetric water content 9. For the second 
and third categories, with the possibilities of changing in the structures of porous media 



and bulk density, the descriptions of these systems may be Eqs. (fTSj) . (J2U) or nonlinear 
fractional diffusion equations described in the text, with the possibilities of using more 
complex diffusivity functions. We note that for the fractional diffusion equation on time the 
value of q = a/2 is less than 1/2, while for the fractional diffusion equation on distance the 
value of q = 1/(1 + a) is greater than 1/2, for < a < 1. These both equations can cover 
the range < q < 1. 



CONCLUSION 



We have investigated the scaling relation (0J) for several types of nonlinear fractional 
diffusion equations. This relation indicates that the nonlinear fractional diffusion equations 
may play important role to the transport of water in soils when the Boltzmann scaling 
is violated. Further, we have shown that the usual nonlinear diffusion equation with the 
diffusion coefficient dependent on the time (|T5j) is closely related to the nonlinear fractional 
diffusion equation (fTTj) . We hope that this work may contribute to reveal some aspects 
related to the fractional diffusion equations and the interplay between fractional and usual 
diffusion equations. 



8 



[1] H. Scher and E. Montroll, Phys. Rev. B 12, 2455 (1975); U. Even et al., Phys. Rev. Lett. 52, 
2164 (1984). 

[2] M. Cates, Phys. Rev. Lett. 53, 926 (1984); H. Schiessel and A. Blumen, Fractals 3, 483 (1995). 
[3] G. Niklasson and C. Granqvist, Phys. Rev. Lett. 56, 256 (1986). 

[4] B. Berkowitz and H. Scher, Phys. Rev. Lett. 79, 4038 (1997); A. Klemm, H. Midler and R. 
Kimmih, Phys. Rev. E 55, 4413 (1997). 

[5] F. Mainardi, Fractinal calculus: Some basic problems in continuum and statistical mechanics, 
in : A. Carpinteri, F. Mainardi (Eds.), Fractals and Fractional Calculus in Continnum Me- 
chanics, Springer- Verlag, NY, 1997, pp. 291-348; T.M. Atanackovic and B. Stankovic, Journal 
of sound and vibration (2004), in press. 

[6] N. Laskin, Physica A 287, 482 (2000). 

[7] D. Hillel, Soil and Water: Physical Principles and Processes, 1971, Academic Press, NY. 
[8] W. Gardner and J.A. Widtsoe, Soil Sci. 11, 215 (1921); Nielsen et al., Soil Sci. Soc. Am. Proc. 
26, 107 (1962); S.I. Rawlins and W.H. Gardner, Soil Sci. Soc. Am. Proc. 27, 507 (1963); H. 
Ferguson and W.R. Gardner, Soil Sci. Soc. Am. Proc. 27, 243 (1963). 
[9] Y. Pachepsky and D. Timlin, J. Hydrology 204, 98 (1998). 
[10] LA. Guerrini and D. Swartzendruber, Soil Sci. Soc. Am. J. 56, 335 (1992). 
[11] R. Metzler and J. Klafter, Phys. Rep. 339, 1 (2000). 

[12] K. B. Oldham and J. Spannier, The fractional calculus (Academic Press, NY) 1974. 

[13] H. B. Callen, Thermodynamics and an introduction to thermostatistics (second edition), 1985, 

John Wiley & Sons. 
[14] M. Bologna, C. Tsallis and P. Grigolini, Phys. Rev. E 62, 2213 (2000). 
[15] J. Klafter, A. Blumen and M. F. Shlesinger 35, 3081 (1987). 

[16] M. C. Jullien, J. Paret and P Tabeling, Phys. Rev. Lett. 82, 2872 (1999); G. Boffetta and 
A. Celani, Phys. Rev. E 60, 6734 (1999); G. Boffetta and I. M. Sokolov, Phys. Rev. Lett. 88, 
094501 (2002). 



9 



